Mesoscopic gap fluctuations in an unconventional superconductor 
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We study mesoscopic disorder fluctuations in an anisotropic gap superconductor, which lead to the spatial 
variations of the local pairing temperature and formation of superconducting islands above the mean-field tran- 
sition. We derive the probability distribution function of the pairing temperatures and superconducting gaps. It 
is shown that above the mean-field transition, a disordered BCS superconductor with an unusual pairing sym- 
metry is described by a network of superconducting islands and metallic regions with a strongly suppressed 
density of states due to superconducting fluctuations. We argue that the phenomena associated with mesoscopic 
disorder fluctuations may also be relevant to the high-temperature superconductors, in particular, to recent STM 
experiments, where gap inhomogeneities have been explicitly observed. It is suggested that the gap fluctuations 
in the pseudogap phase should be directly related to the corresponding fluctuations of the pairing temperature. 
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Understanding the phase diagram and the properties of the 
high-temperature and other unconventional superconductors 
has been among the most complex problems of modern con- 
densed matter physics. Most current theoretical approaches 
to the problem concentrate on strong correlation physics and 
usually assume that the effects of disorder are unimportant. 
However, there exist a number of recent experimental works, 
in particular STM studies of the high-T c cuprateS j 1 ' 2 ' 3 ' 4 ' 5 ' 6 ' 7 
which provide a tentative indication that at least in some ma- 
terials disorder plays an important role in the local formation 
of the superconducting gap. In particular, Gomes et al. 2 have 
studied the local development of the gap as a function of tem- 
perature in Bi2Sr2CaCu20g + 5 above the superconducting tran- 
sition and up to a pseudogap temperature, where the gap in- 
homogeneities cease to exist. An important result of this ex- 
periment is that the real-space gap map observed was static 
and reproducible. This strongly suggests that the inhomoge- 
neous gap formation is unlikely to be a phase-separation or 
superconducting fluctuation effect, but is due to some kind of 
disorder in the system. 

Motivated by these experiments, we theoretically consider a 
disordered superconductor with an unusual pairing symmetry 
(e.g., a (/-wave superconductor) and study mesoscopic vari- 
ations of the local pairing temperature. We point out that 
the existence of the Griffiths-type 8 ' 9 ' 10 phase in the super- 
conducting phase diagram is specific to an anisotropic gap 
superconductor and should not occur in the conventional s- 
wave systems (due to Anderson theorem), unless they are ex- 
tremely dirty or time-reversal symmetry is broken . 11 ' 12 ' 13 An 
important observation is that if the pairing gap is anisotropic, 
the Anderson theorem breaks down and the superconduct- 
ing pairing temperature, T p , is suppressed by disorder even 
if time-reversal symmetry is preserved (here and below we 
make a distinction between the pairing temperature, T p , and 
the superconducting transition temperature, T c , although in 
the framework of the weak-coupling BCS theory they are es- 
sentially the same). The impurities are positioned randomly 
in space and their density is a random variable. Thus, there 
always exist regions where the distribution of impurities is 
such that the local pairing temperature, T p {r), is larger than 
the system-wide average value, (T p ), and the experimental 
temperature, T. 14 These regions form islands with a well- 



defined gap, which exist on the background of a metal, if 
T p (r) > T > (T p ). The width of the "mesoscopic fluctu- 
ation" region certainly depends on the strength of disorder 
and the only parameter, which may enter this dependence, is 
the dimensionless conductance. This defines a narrow win- 
dow where the disorder-induced Griffiths phase co-exists with 
strong superconducting fluctuations. Therefore, the picture of 
impurity induced inhomogeneities in an anisotropic gap BCS 
superconductor is that of superconducting islands and metallic 
regions with strongly suppressed density of states. 

We start with the following Hamiltonian with a built-in l- 
wave pairing (Z > 0): 



<H 



V 2 

2m 



<A(r)-^(r)fc(r) ,(1) 



where U(r) is a disorder potential, A/ is the /-wave interac- 
tion constant, and b(r) corresponds to an "Z-wave Cooper pair" 
b(r) = 2 Xi($)$(k + q/2)<A(p - q/2)e iq r , with $ being the an- 

gle between the direction of the vector k and the x-axis and 
Xi(<P) is the function, which enforces the Z-wave symmetry of 
the gap in the mean-field. 

The first step is to integrate out the fermions and express the 
action in terms of the order parameter Ak = 2 VQt, k')F(k' - 

k,k' 

q/2, k + q/2)e' qr , where according to Eq. ([]]) the interaction 
V(k, k') = —Aixi(<f>)xi(</>') an d F is the standard Gor'kov's 
Green's function. In what follows we will concentrate on the 
spatial dependence of the gap and assume that its symmetry 
in the k-space is preserved: Ak(r) = A(r)^/(k) = Ao/(r)^/(k). 
Using these notations, we arrive at the following free energy 
for the system expressed in terms of the inhomogeneous order 
parameter A(r) 

r[A,C/] = i f A*(ri)A(ri,r2)A(r 2 ) (2) 

1 Jl,2 

+1 f A*(r 1 )A*(r 2 )B(r 1 ,r2,r 3 ,r4)A(r 3 )A(r 4 ), 

4 J 1,2,3,4 

where A(ri,r 2 ) - Xl x 8{t\ - r 2 ) - C(ri, r 2 ) and the Cooperon, 
C, is a random matrix expressed through the Green's func- 
tions (before averaging over disorder) as follows C(ri,r 2 ) = 
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<5C(r 1 ,r 2 )SC(r3,r 4 )> = 



FIG. 1 : (a) The particle-particle bubble for an anisotropic gap su- 
perconductor. This Cooperon diagram contributes to the coefficient 
in the quadratic term of the Ginzburg-Landau expansion ((2). All dis- 
order vertex corrections to the Cooperon vanish.; (b) Pictorial repre- 
sentation for the coefficient in quartic term of the Ginzburg-Landau 
expansion (c) One of the UCF-type dia grams, which contribute 
to the mesoscopic fluctuations of the transition point and mesoscopic 
gap fluctuations. 



r 2 G(e„; ri, r2)G(-e„; ri, r2). As long as we are interested 

in the location of the classical (finite-temperature) phase tran- 
sition, the dynamics of the order parameter and the Cooperon 
are not important. We present the Cooperon as a superposition 
of a local "mean-field" part and a disorder dependent correc- 
tion C = (C) + 6C. The "mean-field" part is diagrammatically 
described by a simple Cooper bubble, without the disorder 
ladder (see Fig. la). Any disorder vertex correction vanishes 
due to the unusual symmetry of the gap. The line where the 
average Ginzburg-Landau coefficient vanishes (A) = deter- 
mines the mean-field transition and leads to the well-known 
Abrikosov-Gor'kov's equation 15 



1 



AnT n T) M2 



(3) 



where Tpo is the pairing temperature without disorder and r 
is the scattering time. This equation implies that the pair- 
breaking effect of the conventional disorder potential in an 
anisotropic gap superconductor [i.e., J d<px(<P) = 0] is iden- 
tical to that of a time-reversal perturbation in an s-wave 
superconductor. 16 We reiterate that Eq. (0 is a result of the 
averaging over disorder in the sample. Below we study meso- 
scopic corrections to this result, which qualitatively can be in- 
terpreted as local changes in the scattering time r(r) in Eq. (0. 

We note that strictly speaking the nonlinear operator in the 
quartic term of Eq. (0 is also random, however its fluctuations 
can be neglected near the transition and its mean field value 
can be used (B). This coefficient is pictorially described by the 
Hikami-box diagram in Fig. lb. A straightforward calculation 
of this diagram gives the following general result: 



(B) 



l6n 2 T 2 



^'"{\ +a ) + ^"[\ +a 



(4) 



where v is the density of states, a - (2jiT p tY 1 , and the over- 
line implies averaging over the Fermi surface. In the clean 
limit, Eq. reproduces the result of Feder and Kallin 17 , 
(B) = 7£(3)v/(8n 2 T 2 ). We note that the dirty limit is not 
reasonable in the context of an anisotropic gap superconduc- 
tor, since it implies that the pairing temperature is suppressed 
to zero and there is no superconductivity. The maximum 
impurity concentration which allows for superconductivity 
(the quantum critical point) is T p uTqpt = y/n ~ 1, where 
y ~ 1.781 is the exponential of the Euler's constant. 

To find the local variations of the transition temperature, 
we consider the following eigenvalue problem for the random 
matrix 6C 



d 2 r'6C(r, r')A(r') = eA(r) 



(5) 



and define the probability distribution function (PDF) of its 
eigenvalues p(s) = (s \e - e[dCVn . The averaging is per- 
formed over the PDF of the random Cooperon matrix, which 



-{6C*k- 



*6C 



where 



we assume Gaussian P |#c] oc exp 
the asterisk implies a convolution over the two spatial vari- 
ables and the operator K corresponds to the correlator of two 
Cooperon operators, which in position representation has the 
form ^r(ri,r 2 ;r 3 ,r4) = (<5C(ri,r 2 )<5C(r3,r4)). This disorder- 
averaged correlator can be calculated using the standard dia- 
grammatic technique (see Fig. lc). These diagrams are topo- 
logically equivalent to the universal conduction fluctuation 
(UCF) diagrams. However there are important differences: 

(i) First, here we are interested in the Cooper channel and 

(ii) Second, we are interested in the local physics, not in a 
long-wavelength behavior of the correlator. 

Using the standard technique, 18 we find the following ex- 
pression for the correlator 



ffi[{r,}] = £(ri-r 4 )S(r 2 -r 3 ) 



An 2 Dv 



(6) 



CO CO 



dt i dt2 



hh(h +t 2 ) : 



exp 



h +h 
4Dt x t 2 



|ri-r 3 |- 



Here the index "1" implies that we consider only one among 
all possible UCF-type diagrams. However, they all contribute 
equally to the correlator of interest and lead to a combinatorial 
factor of c, which is equal to c — 12 in the orthogonal ensem- 
ble and c — 6 in the unitary ensemble (e.g., in the presence of 
a magnetic field). 

The PDF of the local transition temperatures and the cor- 
responding gap amplitudes can be obtained using the optimal 
fluctuation method 19 



<p(e)> oc exp 



(7) 



where the eigenvalue e has the physical meaning of a local 
pairing temperature fluctuation and /(r) is a normalized func- 
tion, which describes the spatial profile and the shape of a sin- 
gle disorder-induced superconducting [if e > (T - (T p ))/T] or 



FIG. 2: (Color online) Plotted are distribution functions of the su- 
perconducting gap A for various temperatures, T. We assumed the 
following parameters: g = 10, E F /T p ~ 40, T p = 93 K. The latter 
choice is motivated by the experimental work 2 on Bi 2 Sr 2 CaCu208 +5 , 
where inhomogeneities have been observed. The temperature behav- 
ior of the PDF following from the mesoscopic fluctuation theory is 
qualitatively similar to that observed in experiment. 2 However, we 
note here that since our approach is based on the BCS theory, the 
theoretical results may not provide a quantitatively accurate descrip- 
tion of the cuprates. 

metallic [if e < (T ~(T p ))/T] puddle. Strictly speaking the lat- 
ter function must be found from a non-linear integral equation 
e/(r) = A f 123 f*( ri )r(r2)K(rur 2 ,r 3 ,r)f(r 3 ) (where A is a 
Lagrange multiplier which appears in the optimal fluctuation 
method; see Ref. [19] for technical details). However, one can 
get a quantitatively reliable description of the PDF by consid- 
ering the puddle function to be a Gaussian of a characteristic 

size i.e., f(r) = (^ 2 ) ex p( - ^)- ^ n principle, one can 
study the distribution of puddle shapes by decomposing the 
function /(r) into spherical harmonics. We do not attempt a 
study of the puddle shapes here, but just point out that "higher- 
orbital momentum puddles" are less probable than spherically 
symmetric ones; the probability of finding a droplet with "mo- 
mentum" m scales as p m oc where po is the probability of 
a spherical puddle. To find the latter we explicitly calculate 

the correlator in Eq. $7$ and find // <£> / |lj / ® /\ = || |r jr, 

where £ is the size of a puddle, I is the mean free path, and 
g = Epr/n is the dimensionless conductance. This leads to 
the following PDF of 7ys ©: 



(pi^Tp)) oc exp -^(y) 9 1 



(T P ) 



(T P ) 



(8) 



We note that by applying a magnetic field, one can cross 
over from the orthogonal to the unitary ensemble and change 
the combinatorial factor in ([8} from c — 12 to c = 6, which 
may be experimentally testable and should manifest itself as 
a diminishing of the random T p (H) or gap distribution width 
exactly by the factor of two. To find the PDF of the supercon- 
ducting gaps, we can just use the Ginzburg-Landau equation 

and set A = ^v(t p - t) /{{T p ){B)), where (B) is given 

by Eq. (0]>. We note that it makes sense to consider a finite- 
size droplet with a well-defined local transition temperature 
or a gap, only if the size of the droplet is much larger than the 



Concentration of dopant atoms x [a.u.] 

FIG. 3: (Color online) This figure is to illustrate the qualitative 
discussion in the text about the possibility of a dome-shaped dop- 
ing dependence of the pairing temperature within the Abrokosov- 
Gor'kov theory. The figure shows a typical dependence of T p on 
"doping," x, in the toy model of the Abrikosov-Gor'kov theory with 
both T pa (x) and t(x) being dependent on the same doping param- 
eter. The graphs are solutions of the following equation for 
ln^x 1- "/^) = 1^(1/2 + x5// p ) - t/f(l/2), where t p is a dimension- 
less pairing temperature and S characterizes the strength of disorder. 
The graphs correspond to a = 1/2 and 6 = 0. 1 (red line), 0.15 (blue 
line), and 0.2 (green line). 



coherence length £ » £a . The opposite limit corresponds to 
the case of a mesoscopic superconducting nanograin in which 
the notion of the gap is not well-defined (see Ref. 12011 for a 
review). Therefore, the smallest possible size of the puddle 
with a well defined gap Ao is £ m ;„ ~ v^/Aq, which implies that 
the dimensionless ratio % m i n /l ~ (T p t) . The latter parameter 
is of order one and thus the PDF of the gaps takes the form 
(see also Fig. [2]) 



P[A] 



2g^bA 
^(T„) : 



■ exp 



bA 2 T 
+ 



(T P ) 



(Tp) 2 



(T P ) 



(9) 



where in the of-wave case, b ~ 
-[%r'{\+») + 2>r{\+<y)\l02n 2 \ a 

{E F l{T p )){2n 2 gY\ and B ~ 1. Note that in Eq. © we 
have omitted a term proportional to oc 6(A), which describes 
normal regions. The physical picture, which emerges for such 
a disordered superconductor is that right above the mean-field 
transition temperature there exist rare superconducting islands 
separated by "normal regions" (which are still very close 
to the local transition temperature). We note that since the 
parameter gT p o/Ep is at best of order one, the "mesoscopic 
Griffiths phase" overlaps with the Ginzburg region of strong 
superconducting fluctuations.— This leads to the conclusion 
that the Griffiths phase is a mixture of superconducting 
islands and metallic regions with strongly suppressed density 
of states. 

Even though our quantitative description directly applies 
only to a weakly coupled BCS superconductor, we believe 
that some aspects of the theory are relevant to the cuprates as 
well ( the importance of disorder effects for the cuprates have 
been discussed previously, see, e.g. [ 22 23 24 25. 26U27D). But 



4 



first, we make the following curious observation: In a high- 
T c superconductor the major source of disorder is presum- 
ably the dopant atoms. But this access oxygen is also the 
source of the carriers, which lead to superconductivity in 
the first place. Thus, the "clean" pairing temperature in the 
Abrikosov-Gor'kov's formula (01 should explicitly depend on 
the doping level, T p o(x) (e.g. through a BCS-like high-energy 
cut-off or possibly via a locally modulated electron pairing 
interaction. 25 ) The scattering time depends on x too and con- 
tributes to superconductivity suppression. Thus there are two 
competing effects of the dopants: They both enhance and sup- 
press superconducting properties. An interesting result, which 
follows from the simple Abrikosov-Gor'kov's equation (0, 
is that even if T p q(x) is monotonically increasing with x, but 
T p o(x)t(x) oc x~ a with a > 0, then the actual pairing tempera- 
ture doping dependence has a dome-shaped form. An example 
of such a dependence for a = 1/2 is plotted in Fig. [3] 

The mesoscopic disorder fluctuations in the Abrikosov- 
Gor'kov theory may not correspond to the modulations of the 
transition temperature in a high-r c superconductor, but should 
be related to the modulations of the pseudogap temperature, 
T t , which is believed to be the onset of Cooper pairing (i.e., 
to T p , but not T c ). The "pseudogap" region T c < T < T p pre- 
sumably represents the regime of strong phase fluctuations 28 
and the superconducting transition is expected to be that of 
XT-type. In the latter scenario, the transition temperature is 
proportional to the superfluid density, which in turn is directly 
related to the local value of the gap. The mean field gap is 



determined by the deviation (T p - T) from the local pairing 
temperature (e.g., via a non-linear BCS-like self-consistency 
equation 6S/6A - 0). If T p (r) fluctuates due to disorder, 
so does A(r) and, quite generally, the local gap should "fol- 
low" local T p . In fact, such a correlation has been observed 
in experiment. 2 In the model of uncorrected short range dis- 
order, the only possible result for the corresponding PDF is 
Eq. (0, with A centered in the vicinity of the mean-field gap 
at the given temperature P[A] oc exp{-/3g 2 [1 - A/<A(7/))] 2 }. 
In a clean material, the corresponding regime of mesoscopic 
fluctuations is very narrow unless there are other types of dis- 
order effects, such as "structural disorder" (e.g., extended de- 
fects, warping of the 2D planes, etc.), which can be modelled 
as a random diffusion coefficient 13 (D(r)D(O)) = (D) 2 a 2 5{r). 
These phenomena lead to qualitatively the same effect of ran- 
dom T p and A, but, may occur at the length-scales much 
larger than the mean-free path and become important in a 
much wider range of parameters [the width of the distribu- 
tion is determined by (a/t)g instead of g ]. We also note 
that if indeed the local gap fluctuations observed in experi- 
ment 2 are related to mesoscopic disorder effects, they should 
be correlated with the pinning properties, 29 i.e., the width of 
the gap distribution should be proportional to the critical cur- 
rent j c in the collective pinning regime 13 , ^(1 - A/(A)) 2 ^ oc 
(jcl jco)(H/H C 2) (where j c o is the critical current in zero field). 

The author acknowledges the Aspen Center for Physics for 
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